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We consider the dynamics of a causal bulk viscous cosmological fluid filled flat constantly decel- 
erating Bianchi type I space-time. The matter component of the Universe is assumed to satisfy a 
linear barotropic equation of state and the state equation of the small temperature Boltzmann gas. 
The resulting cosmological models satisfy the condition of smallness of the viscous stress. The time 
evolution of the relaxation time, temperature, bulk viscosity coefficient and comoving entropy of the 

- - - dissipative cosmological fluid is also obtained. 
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. I. INTRODUCTION 

» , 

o . 

' In a series of recent papers - |3|, a viscous Bianchi type I Universe with metric of Kasner form has been 
investigated. The cosmic fluid was supposed to be endowed with a bulk viscosity ^ and a shear viscosity rj, with 
both viscosity coefficients independents of position. The corresponding energy-momentum tensor is Tit = pUiUk + 
' [p ~ ^9) hik — 2rjaik , with Ui the four -velocity, htk = gik + UiUk the projection tensor , = u\ the scalar expansion and 
<^ik — h {ui-ih''), + Uk;ih\^ ~ \hikQ the shear tensor. This energy-momentum tensor is of Eckart-Landau (non-causal) 
^ \ type [Q, The entropy current four- vector is given, in these type of models and in the absence of the heat flux, 
' by 5*^ = nkB<7u^ , where ks is the Boltzmann constant, n is the baryon number density and a the nondimensional 
' entropy per baryon. But as has been recently shown in |^ , if the thermodynamics together with the dominant energy 
, conditions are kept, the viscous cosmological fluid model in a Kasner type geometry cannot describe the growth of 
^—1 ' the entropy in the Universe. 

T-H , In many cosmological and astrophysical situations an idealised fluid model of the matter is inappropriate. Such 
■ possible situations are the relativistic transport of photons, mixtures of cosmic elementary particles, the evolution of 
^ . cosmic strings due to their interaction with each other and with the surrounding matter, the classical description of 
O^' the (quantum) particle production phase, interaction between matter and radiation, quark and gluon plasma viscosity, 
^' ^ ^ different components of dark matter etc. From a physical point of view the inclusion of dissipative terms in the 
bJ[)' energy momentum tensor of the cosmological fluid seems to be the best-motivated generalisation of the matter term 
^ ' of the gravitational field equations. 
• ^ , The theories of dissipation in Eckart-Landau formulation Q , |^ , who made the first attempt at creating a relativistic 
' theory of viscosity, are now known to be pathological in several respects. Regardless of the choice of equation of state, 
^ , all equilibrium states in these theories are unstable. In addition, as shown by Israel 0, signals may be propagated 
■ - - ' through the fluid at velocities exceeding the speed of light. These problems arise due to the first order nature 
of the theory since it considers only first-order deviations from the equilibrium, leading to parabolic differential 
equations, hence to infinite speeds of propagation for heat flow and viscosity, in contradiction with the principle of 
causality. While such paradoxes appear particularly glaring in relativistic theory, infinite propagation speeds already 
constitutes a difficulty at the classical level, since one does not expect thermal disturbances to be carried faster than 
some (suitably defined) mean molecular speed. Conventional theory is thus applicable only to phenomena which are 
"quasi-stationary" i.e. slowly varying on space and time scales characterised by mean free path and mean collision 
time 0. This is inadequate for many phenomena in high-energy astrophysics and relativistic cosmology involving 
steep gradients or rapid variations. These deficiencies can be traced to the fact that the conventional theories (both 
classical and relativistic) make overly restrictive hypothesis concerning the relation between the fluxes and densities 
of entropy, energy and particle number. 

A relativistic second-order theory was found by Israel Q and developed by Israel and Stewart Q into what is called 
'transient' or 'extended' irreversible thermodynamics. In this model deviations from equilibrium (bulk stress, heat 
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flow and shear stress) are treated as independent dynamical variables leading to a total of 14 dynamical fluid variables 
to be determined. The solutions of the full causal theory are well behaved for all times. Hence the advantages of the 
causal theories are the foUowings 1) for stable fluid configurations the dissipative signals propagate causally 2) 
unlike Eckart-type's theories, there is no generic short-wavelength secular instability in causal theories and 3) even 
for rotating fluids, the perturbations have a well-posed initial value problem. Therefore, the best currently available 
theory for analysing dissipative processes in the Universe is the full Israel-Stewart causal thermodynamics. 

Due to the complicated nonlinear character of the evolution equations, very few exact cosmological solutions of the 
gravitational field equations are known in the framework of the full causal theory. For a homogeneous Universe filled 
with a full causal viscous fiuid source obeying the relation £, ^ , exact general solutions of the field equations have 
been obtained in and - . In this case the evolution of the bulk viscous cosmological model can be reduced 
to a Painleve-Ince type differential equation. It has also been proposed that causal bulk viscous thermodynamics can 
model on a phenomcnological level matter creation in the early Universe , ||l^ . 

Because of technical reasons, most investigations of dissipative causal cosmological models have assumed FRW 
symmetry (i.e. homogeneity and isotropy) or small perturbations around it |^5|. The Einstein field equations for 
homogeneous models with dissipative fluids can be decoupled and therefore are reduced to an autonomous system of 
ordinary differential equations, which can be analysed qualitatively [ p^ , [T^ . A qualitative analysis of the Bianchi type 
I cosmological model has been performed, under the assumption of the plane symmetry, by Romano and Pavon [[l8[ . 
Homogeneous and isotropic cosmological fluids with bulk viscous dissipation have been analyzed, in the framework of 
the truncated Israel-Stewart theory, with the use of dynamical systems techniques, by Di Frisco, Herrera and Ibanez 
p9[ . They found that almost all solutions inflate and only few of them can be considered physical, since the dominant 
energy condition is not satisfied. 

The matter density parameter Q and the deceleration parameter q are the most important observational parameters 
in cosmology. The values of the deceleration parameter separates decelerating (q > ) from accelerating { q < 0) 
periods in the evolution of the Universe. Determining the deceleration parameter from the count magnitude relation 
for galaxies is a difficult task due to evolutionary effects. The tightest limits on the present value qo of the deceleration 
parameter obtained from observations are —1.27 < f/o < 2 po|| . Studies of galaxy counts from redshift surveys provide 
a value of go = O-lj with an upper limit of go < 0-75 pO[ |. Recent observations show that the deceleration parameter 
of the Universe q is in the range — 1 < g < 0, and the present-day Universe undergoes an accelerated expansionary 
evolution |2^]- But of course these results do not exclude the existence of a decelerating phase in the early history 
of our Universe. 

Cosmological models with constant deceleration parameter have been intensively investigated in the physical litera- 
ture. It has been shown by Berman and Gomide ||2J] that all the phases of the Universe, i.e., radiation and pressure-free 
phase, may be considered as particular cases of the q ~ constant type. Incidentally, most of the well-known perfect 
fiuid models in general relativity and Brans-Dicke theory are models with constant deceleration parameter. 

It is the purpose of the present paper to consider causal bulk viscous cosmological models in a fiat Bianchi type I 
space-time. We assume that during the cosmic evolution the deceleration parameter is constant and positive. Thus 
we extend to the anisotropic case the previous investigation of isotropic fiat and homogeneous constantly decelerating 
bulk viscous cosmological models |2^. This model satisfies the condition of smallness of the viscous stress. We 
also analyse the evolution of the physical parameters of the dissipative cosmological fluid such as relaxation time, 
temperature, bulk viscosity coefficient and comoving entropy. For the equilibrium pressure we shall adopt a linear 
barotropic equation of state and the state equation of the small temperature Boltzmann gas. 

The present paper is organised as follows. The basic equations describing the dynamics of a decelerating causal bulk 
viscous cosmological model are obtained in Section II. The case of a viscous cosmological fluid obeying a barotropic 
equation of state is considered in Section HI. Boltzmann gas flUed viscous Bianchi type I Universes are discussed in 
Section IV. In Section V we discuss and conclude our results. 



II. THERMODYNAMICS, FIELD EQUATIONS AND CONSEQUENCES 

The energy momentum tensor of a relativistic fluid, with bulk viscosity as the only dissipative phenomenon is 

Tl" = ip + p + n)u,u'' -ip + U)5f. (1) 

In Eq. (Q) p is the energy density, ,i = 0, 1, 2, 3 , is the four-velocity, ( Uiu'^ — 1), p is the equilibrium pressure 
and n is the bulk viscous pressure. The particle flow vector is given by — nu^, where n > is the particle number 
density. N\ = leads to the particle number conservation equation h + ZnH ~ 0, with H = the Hubble factor 
(a semicolon denotes the covariant derivative with respect to the metric). The conservation law T*^. — implies 
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p + 3F (p + p) = -3HII. (2) 

In the framework of causal thermodynamics and hmiting ourselves to second-order deviations from equilibrium, the 
entropy flow vector takes the form | |2^ : 

^ sN^ + = sN^ - —u\ (3) 

where s is the entropy per particle, r is the relaxation time, T is the temperature and ^ is the coefficient of bulk 
viscosity. 

The causal evolution equation for the bulk viscous pressure 11 is given by 

rn + n = -i^H - irn I 3i/ + 3 - i - . (4) 

Equation arises as the simplest way (linear in 11) to satisfy the H theorem (i.e., for the entropy production 
to be non-negative, 5!^ = |y > 0) 0], |26|. The Israel-Stewart theory js) is derived under the assumption that the 




thermodynamic state of the fluid is close to equilibrium, which means that the non-equilibrium bulk viscous pressure 
should be small when compared to the local equilibrium pressure, that is |n| < p. If this condition is violated then one 
is effectively assuming that the linear theory holds also in the non-linear regime far from equilibrium but for a fluid 
description of the matter, the condition ought to be satisfied. Due to the non-thermalizing character of dissipative 
expansion which require that the rate of interactions maintaining equilibrium remains lower than the expansion rate 
we have the following consistency condition on causal viscous cosmologies: - < H. 

We assume equations of state in the general form p — p {n, T) ,p — p (n, T), according to which the particle number 
density n and the temperature T are the basic thermodynamical variables. In this case one finds the following 
evolution law for the temperature of a causal bulk viscous cosmological fluid 

where we denoted b = (dp/dT) / (dp/dT) ,a^p/ {Tdp/dT) and i = {p + p) I p. 

The growth of the total comoving entropy over a proper time interval (toj^) is given by p6t 

S(t)-S{to)^-^J^'^dt, (6) 

where V is the proper comoving volume element. 

In the case of a flat Bianchi type I space-time with a line element 

ds'^ = dt'^ - al{t)dx^ - al{t)dy^ - al(t)dz^ , (7) 

the gravitational field equations take the form 

iH + Hl + Hl + Hl^-]^{p + ip + m), (8) 



^^(^^ff,:) = ^(p-p-n),* = l,2,3, (9) 



where we denoted V — aia2a3 , Hi = = 1,2,3 and H = ^ J2'i=i ~ sv 
By adding Eqs. (|) we obtain 
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leading to 



H + SH^^lip-p-U), (10) 



H, = H + ^,t = l,2,3, (11) 
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where Ki,i = 1,2,3 are constants of integration satisfying the consistency condition J2'i=i — ^■ 
The energy density p of of the bulk viscous cosmological fluid is given by 

P-^H^i^-^). (12) 

where K'^ = Y^l^^Kf. The weak energy condition p > imposes the constraint < QH^V^ on the integration 
constants Ki, i = 1,2, 3. 

As a first example of a Bianchi type I space-time we consider the Kasner type Universe for which the components of 
the metric tensor are given by ai{t) =tP%i= 1,2,3 and pi = const.. We denote M = J2^=i Pi- Hence we have V = 
, Hi = = 1,2,3 and H = The coefficients pi,i = 1, 2, 3 must satisfy the consistency conditions, which follows 
from ^ {yHi) = H + 3H^ , i = 1,2,3 and are given by 3pi = Pi + P2 + Ps^i = 1,2, 3. These equations represent a 
homogeneous system of algebraic equations. Since the determinant of the system is zero, it admits a unique non-trivial 
solution given by pi = p^,p2 = Pa- Consequently, it follows that there are no anisotropic Kasner- type causal bulk 
viscous cosmological fluid flUed space-times. On the other hand this result is generally valid for all cosmological fluids 
with isotropic pressure distribution. 

A useful observational quantity in cosmology is the deceleration parameter 

The positive sign of q corresponds to standard decelerating models whereas the negative sign indicates inflation. 

For a barotropic (p = (7—1)^,7 = const, and 1 < 7 < 2) cosmological perfect fluid fllled flat isotropic FRW 
type space -time the deceleration parameter is a constant for all 7 . In this case the gravitational field equations are 
3H^ = p, 2H + 3H^ = — (7 — 1) p, leading to q = ^ — 1, respectively. The deceleration parameter is constant in 
both radiation (7 = |) and matter (7 = 1) dominated phases. For a barotropic perfect fluid fllled Bianchi type I 
geometry using Eq. ( pJ3[ ) the deceleration parameter can be expressed as 

,-2 (2-7)p _37 ^ 9(2^7) i^^ .,4. 
^" 2H^ ~ 2 ^ 4 V2- ^^^> 

For 7 = 2 (stiff or Zeldovich fluid) q is constant for all times. For a rapidly expanding Universe with V >> K , it 
is a very good approximation to also consider the deceleration parameter as a constant (g is an exact constant for 
V = const. ). 

Since the bulk viscous pressure represents only a small correction to the thermodynamical pressure, it is a reasonable 
assumption that the inclusion of viscous terms in the energy -momentum tensor does not change fundamentally the 
dynamics of the cosmic evolution. Therefore in the foUowings we restrict our analysis to cosmological evolutions 
characterised by the constancy of the deceleration parameter, 

q = qo = constant > 0. (15) 

With this assumption the Hubble factor H immediately follows from Eq. ( p^ ) and is given hy H ^ where 
Ho = = constant. Hence the volume factor is V = Vot^^° , where Vq > is a constant of integration. With the 
Hubble factor fixed by Eq. ( p^ ) the time evolution of the scale factors is given by 

a,{t) ^ aMt^'e^ou-sHoji"''"" , i l, 2, 3, (16) 

where a^o, i = 1,2,3 are constants of integration. 

From the field equations we obtain the energy density of the cosmological fluid in the form 

Po A C 



Pit)-f[^-jm^)' (17) 

where pq = 3Hq and C = ^^^^ > 0. The sum of the thermodynamic and viscous pressure is 

Hoi2-3Ho) / 1_\ 

t^ { 2Hoi2^3Ho)VoH^"o-2j- 



4 



For 6-ffo — 2 > 0, or, equivalently, qq < 2, the energy density p of the cosmological fluid is positive only for time 



intervals t > to ~ { Q^j-iy'i 



We also introduce the anisotropy parameter defined by 

The variation of the particle number density of the bulk viscous cosmological fluid follows from n — ^ — -^§^, 
where Uq > and no > are constants of integration. 

The general solution of the gravitational field equations for a Bianchi type I space-time with causal bulk viscous 
fluid can be generally expressed in an exact form for an arbitrary time dependent deceleration parameter q = q{t). In 
the Appendix we present the general representation of the solution in terms of the deceleration parameter. 



III. MODELS WITH A LINEAR BAROTROPIC EQUATION OF STATE 



In order to completely solve the gravitational field and the causal evolution Eqs. (||) and (^), together with the 
condition of the constancy of the deceleration parameter for an anisotropic Bianchi type I geometry, we must close 
the system by an equation of state for p and specify the thermodynamic coefficients ^ and t. We assume first that 
the thermodynamic pressure of the bulk viscous cosmological fiuid obeys a linear barotropic equation of state 



M^) = (7-l)p(t) = ^^^(l- 



C 



(20) 



where 7 (the adiabatic index) is a constant and 1 < 7 < 2. Hence the field equations (|T^) and (^7|) can be immediately 
solved to give: 



Ho 

i2 



(2 - 7) gi 

teHo-2 



(21) 



where IIq = Hq (2 — 87^0) and Ci 



73T 



Since for a viscous type dissipative theory it is expected that the bulk viscous pressure is negative, 11 < 0, it follows 



that 2 — 3^Ho < leading to go < 



37-2 



For 7 = I (radiation gas) we obtain go < 1- 



The causal evolution Eq. (^ can be formally integrated to give 



T(i) = 7^oi^exp 



tB 



-)dt 



(22) 



with Tqi a constant of integration. 



Following Belinskii, Nikomarov and Khalatnikov |27|, we suppose that the relation t ^ p holds in order to guarantee 



that the propagation velocity of bulk viscous perturbations, i.e., the non-adiabatic contribution to the speed of sound 

1/2 

in a dissipative fluid without heat flux or shear viscosity, which is expected to be of the order, Cb = [^/ {p + p) t] ' ~ 

[C/pt]^ does not exceed the speed of light. With this choice for the relaxation time r the causal evolution Eq. (Q) 
gives the following evolution law for the temperature: 



T{t) = Tot' 1 



C 



(2-7)^1 



2 f- 



■dt 



(23) 



where we denoted Tq = Tqi^ , I ^ + 3i/n - 2 and to = 2 



po ' - Uo ■ - — ^ - ■ (2-7)no' 

An analysis of the relativistic kinetic equation for some simple cases given by Murphy [g8| , Belinskii and Khalatnikov 
p9[ and Belinskii, Nikomarov and Khalatnikov has shown that in the asymptotic regions of small and large values 
of the energy density, the viscosity coefficients can be approximated by power functions of the energy density with 
definite requirements on the exponents of these functions. For small values of the energy density it is reasonable to 
consider large exponents, equal in the extreme case to one. For large p the power of the bulk viscosity coefficient 
should be considered smaller (or equal) to 1/2. Hence, we shall assume that the bulk viscosity coefficient obeys the 
simple phenomenological law 
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(24) 



with a and s constants. Consequently, using Eq. (|2j) the evolution of the temperature in a bulk viscous cosmological 
fluid filled constantly decelerating Bianchi type I Universe is given by 



T{t) = Tot' 1 - 



C 



l-(2-7) 



(25) 



where F{t) ~ exp 



,2(3-1) 



dt 



The evolution of the temperature for time intervals t > C^^o-^ is presented, for a radiation dominated cosmological 
fluid (7 I) in Fig.l. 




1.2 



1 . 3 



1 . 4 



FIG. 1. Behaviour of the temperature T{t) of the barotropic causal bulk viscous cosmological fluid for different values of the 
deceleration parameter: qo ~ j (solid curve), go = ^ (dotted curve) and qo ~ j (dashed curve). 



For all times and for values of go < 2 the temperature is a monotonically decreasing function of time. 

The behaviour of the comoving entropy is governed by Eq. (^ . Since the bulk viscous pressure 11 < , from Eq. (|^) 
it follows that |f > 0. The variation of the entropy as a function of time for different values of the deceleration 
parameter is presented in Fig. 2. 




6 



The entropy is a monotonically increasing function of time for deceleration parameters in the assumed range. Hence 
due to bulk viscous type dissipative stresses a large amount of entropy can be produced during the very early evolution 
of our universe. 

Finally, we shall check the condition of the smallness of the non-equilibrium bulk viscous pressure Tl{t) with respect 
to the local equilibrium pressure, p(t), meaning that the thermodynamical state of the fluid is close to equilibrium. 
From Eqs. (fg) and (|l]) we obtain 



U{t) 




Ho 


Pit) 




(7 - 1) Po 



1 



c 



1 - 



3(2 - 7)Cgo 1 



< 1. 



(26) 



The condition ( |2^ ) is satisfied for qo < 2 (in the small time limit) and qq < ^''^ ^ (large time limit). For 7 = | we 
obtain go < §• Therefore, the cosmological solutions obtained above are thermodynamically consistent. 



IV. MATTER DOMINATED CAUSAL BULK VISCOUS BIANCHI TYPE I COSMOLOGICAL MODELS 

In the previous Section we have considered that the bulk viscous cosmological fluid filling the early universe satisfies 
a linear barotropic equation of state, thus having a radiation-like behaviour. We consider now the case of a bulk 
viscous matter dominated cosmological fiuid consisting from a relativistic classical monatomic Boltzmann gas (i.e. 
a dilute monatomic gas with high collision rate). By denoting (3 = the ideal-gas law for the Boltzmann gas 

take the form: p — nmf3^ , nm — Aoe-°'/3-^K2 (/?) and p ^ Aq [P^^Ki {(3) + 3/3^2^2 (/?)] , with a the relativistic 
chemical potential, the constant Aq = m^g /2TTh^ , g is the spin weight of the fluid particles and Ki,i = 1,2,3 are 
modifled Bessel (Hankel) functions. For a Boltzmann gas the relaxation time and the bulk viscosity coefficient can 
be expressed as t = l/ja^ and ^ = mcrj'^il'^ /jq, with -q = {(3) / K2 {(3), ft — 3-f[l + l/(3r]] — 5 and 7 is a solution 
to [7/(7 - 1)] = (3^ (1 + {bril[3) - if) @. 70 is the colhsion integral. 

The equations governing a flat Robertson- Walker cosmological model containing a dissipative Boltzmann gas have 
been integrated numerically within the framework of the full causal dissipative thermodynamics by Hiscock and 
Salmonson In this case the evolution of the Universe is non-infiationary and the pressure is not allowed to become 
negative. Inhomogeneous Lemaitre-Tolman-Bondi type hydrodynamic models for a cosmological matter fluid with 
causal viscous stresses but no heat conduction satisfying the above equations of state have recently been considered 
indl-ll. 

The equation of state of the Boltzmann gas in the form given above are not especially enlightening. Hence it is more 
useful to evaluate the state equation of the relativistic Boltzmann gas in two regimes: in the high temperature limit 
corresponding to >> 1 and in the small temperature limit for which the physical parameters of the gas satisfies 
the condition ^ <<1. 

In the small temperature limit the equation of state of the Boltzmann gas is 

3 T 

p = mn+ -nT,P = nT,— « 1. (27) 
2 m 

For a matter fluid satisfying equations of state of the form (p^ ) we have 7 = 1, a = (2m/3T) >> l,and b — (2/3). 
We assume that the gas evolves in a flat anisotropic Bianchi type I space-time and the cosmological evolution is 
characterised by a constant positive deceleration parameter, q = qo = constant > 0. With this assumption the energy 
density of the low temperature Boltzmann gas is given again by Eq. (p^. Then from the gravitational field equations 
(|)-(|) and with the use of Eqs. (0) and @) we obtain: 

Pit) ^ ^15. 1 - 77^77:^ - „ ' (28) 



3 t'^Ho-'^ po t^H,,- 

^^^^ " 3no i2-3Ho ' Ipi^ ~ ^3^0-2 ) ' (29) 



Hp (2 - Sffo) 

i2 



CHp 1 2 mrio 



2-5Hot^fi'>-2 3 Ho i2 ~ 5Ho) t^H"-^ 



(30) 



Equations (|2^)-(|3^) are physically meaningful only for time intervals satisfying the condition 1 — ^6ul-i 
i-iUo-i > 0, thus leading to positively defined thermodynamic pressure and temperature. 
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The requirement that the temperature be a decreasing function of time leads to the condition 2 — 3Hq > or to 
the restriction 90 > 5 imposed on the deceleration parameter. The condition of the negativity of the bulk viscous 
pressure implies go < §• 

The causal evolution equation (^) for the bulk viscous pressure can be written in the form 

2 r„ , /i3Hon2- 



d 

di 



n \t 



2 d , 

T dt 



T 



leading to the following expression for f = ^: 



/ = 



where C2 is a constant of integration. 



exp(2/f) 



T 



C2 - 6H0 J 



exp (2 / 



dt 



(31) 



(32) 



If the condition 



"r2 n 



>> 



2? 



or, equivalently, |^| >> l^lis fulfilled, we can neglect in Eq. (Is^) the term -^f 



In this case we obtain for |- the following approximate form: 



^ _ t3Hon2 
T 



T 



C2 - 6H0 



^3ffo-in 



T 



-dt 



(33) 



The variation of the bulk viscosity coefficient-relaxation time ratio is represented in Fig. 3. 
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FIG. 3. Time evolution of the bulk viscosity coefficient-relaxation time ratio f{t) — | for the low temperature Boltzmann 
gas filled Bianchi type I space-time for different values of the deceleration parameter: qo = 0.6 (solid curve), qo — 0.8 (dotted 
curve) and go = 1.1 (dashed curve). 

Another thermodynamic model can be obtained if we consider that the relaxation time is related to the energy 
density of the low temperature Boltzmann gas by means of the relation ^ ~ p. 
Then the causal evolution Eq. (^) leads to the general relation 

(34) 



18g3 



A 

dt 



In 



V PT ) 



The behaviour of the relaxation time t is presented, for different values of the deceleration parameter, in Fig. 4. 
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FIG. 4. Dynamics of the relaxation time T{t) in the low temperature Boltzmann gas filled Bianchi type I Universe for 
different values of the deceleration parameter: qo = 1.05 (solid curve), go — 1.15 (dotted curve) and qo = 1.25 (dashed curve). 

In the large time limit, as expected, the relaxation time tends to zero. 

The variation of the comoving entropy of a monatomic gas in a constantly decelerating Universe can be obtained 
from Eq. (|^). In Fig. 5 we present the time evolution of the entropy for the low temperature Boltzmann gas filling an 
anisotropic constantly decelerating Bianchi type I Universe. 
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FIG. 5. Variation of the entropy S{t) of the low temperature Boltzmann gas filled Bianchi type I Universe for different values 
of the deceleration parameter: go = 0.6 (solid curve), go — 0.8 (dotted curve) and go — 1.1 (dashed curve). 



mt) 



< 1. In the 



The condition of the smallness of the bulk viscous pressure can be studied from the relation 

limit of large and small times we obtain the restrictions go < f and qq < 2, respectively. Hence the thermodynamic 
consistency of the low temperature dissipative Boltzmann gas filled Bianchi type I Universe is realised for values of 
the deceleration parameter in the range qo G (^i §)• 

The Boltzmann gas is not a physically realistic model of the content of the Universe at all temperatures and densities 
(among other problems it ignores the creation of particle-antiparticle pairs at relativistic temperatures). But it is 
(together with the radiation fiuid) the best physically motivated and analysed cosmological fluid model. Even the low 
temperature approach is still reasonable at high temperatures (for example, for the electron m w 10^ K and the low 
temperature limit equations of state should be very accurate up to about 10^ K). 
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V. DISCUSSIONS AND FINAL REMARKS 



The study of the evolution and dynamics of the early Universe is usually performed by looking for the consequences 
of the applications of the physical laws to the cosmological environment. In the present paper we have used the 
inverse approach, by investigating the effects of the cosmology on the laws of the physics. At least for small (on 
a cosmological scale) intervals of time the deceleration parameter can be considered a constant or a slowly varying 
parameter. The conditions of the constancy and of the non-negativity of q impose strong constraints on the evolution 
of the cosmological models. For causal bulk viscous cosmological models these restrictions lead, once we adopt an 
equation of state for the dense matter and via the gravitational field equations, to explicit expressions for the energy 
density, thermodynamic pressure and bulk viscous pressure. In the present paper we have considered two distinct 
models, depending on the equation of state of the cosmological fluid and on the relation. In all these cases explicit 
expressions for the time dependence of the thermodynamic parameters (temperature, bulk viscosity coefficient and 
relaxation time) can be obtained. In these physical models the thermodynamic consistency of the results requires 
values of q so that q < | , which is consistent with the range of values of the deceleration parameter obtained from 
observations [ pO| . 

In the large time limit, t ^ oo and if 1 — 3i?o < or, equivalently, go < 2 , from Eqs. (|l6|) it follows that 
ai{t) t^°,i = 1,2,3 and the mean anisotropy parameter yl ^ 0. Hence due to the presence of dissipative stresses 
a Bianchi type I universe ends in an isotropic phase. The transition from the anisotropic to the isotropic state 
is associated to a large increase in the entropy. The increase of the comoving entropy of the cosmological fluid is 
independent, in the causal thermodynamical approach, of the equation of state of the matter and seems to be a general 
feature of the model. Hence this approach can lead to a better understanding of the entropy generation mechanisms 
in the very early Universe, where viscous type physical processes have probably played an important role even when 
the Universe was about 1000s old. 

Therefore the constantly decelerating causal bulk viscous cosmological fluid model in an anisotropic geometry leads 
to a self-consistent thermodynamic description which maybe could describe a well-determined period of the evolution 
of our Universe. 
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APPENDIX 

In this Appendix we present the general solution of the gravitational field equations for a causal bulk viscous 
cosmological fluid filled Bianchi type I Universe in the case of a known functional form of the deceleration parameter 
q — q(t). From the definition of q we obtain first 



H 



qi+t+ I q{t)dt 



(35) 



with qi a constant of integration (in the main text we have obtained all the results with the particular choice qi — 0). 
With the use of 3HV — V the volume scale factor can be represented as 



V = Vo exp 



dt 



qi+t + J q{t)dt 



(36) 



with Vq a constant of integration. The scale factors are 



fli = exp 



dt 



K, 



qi+t + J q{t)dt Vo 



dt 



exp ( 3 f ^ 

I qi+t+ J q(t)dt 



1,2,3, 



(37) 



with a^o, i — 1, 2, 3 arbitrary constants. 

From Eqs.(||) and ( p^ ) we obtain the energy density of the bulk viscous fluid in the general form 

,^3/7^-4=377^(1-^ 



(38) 



We assume that the cosmological fluid obeys the linear barotropic equation of state p = (7 — 1) p. Then the general 
form of the bulk viscous pressure follows from the gravitational field equations and is given by 



n = p-p- 2 



[h + 37f2^ 



(39) 



The condition of the negativity of the bulk viscous pressure is satisfied only if the deceleration q[t) and the mean 
anisotropy A{t) parameters obey the inequality 



2[l + q(t)]-37<-(2-7)A(t). 



(40) 



From Eq. ( P2|) and with the use of the evolution laws r = ap^ ^ and ^ = ap" for the relaxation time and bulk 
viscosity coefficient respectively, we obtain the evolution law of the temperature in the form 



T{t) = Tf 



(2-7)p-2 7r + 377 



01- 







1=/ 









(2 -7)p- 2 77 + 3772 




dt 



(41) 



The condition of the smallness of the non-equilibrium bulk viscous pressure with respect to the local equilibrium 
pressure can be checked from the relation 



n 




p 





2V^ 



- 377^2 - 27r 



3(7- 1)772 (1 -^; 



2 ^v'-y r 3772^ dt 



(7-l)(l-^ 



(42) 
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